Electronic properties of the degenerate Hubbard Model 
A dynamical mean field approach 
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We have investigated electronic properties of the degenerate multi-orbital Hubbard model, in the 
limit of large spatial dimension. A new local model, including a doubly degenerate strongly corre- 
lated site has been introduced and solved in the framework of the non-crossing approximation (NCA). 
Mott-Hubbard transitions have been examined in details, including the calculation of Coulomb re- 
pulsion critical values and electronic densities of states for any regime of parameters. 



An increasing interest for strongly correlated elec- 
tronic systems has been induced by the discovery of high 
temperature superconductors ||l[ and giant magneto- 
resistance [EJ . Despite a large amount of papers on this 
fascinating field, a complete physical understanding of 
strongly correlated electronic systems remains a major 
problem of condensed matter theory. Even the simplest 
Hamiltonian describing such systems, given by the Hub- 
bard model, is still an unsolved challenging problem. 

The seminal work of Metzner and VoUhardt g has 
shown the importance of large spatial dimensions D. 
Taking the limit D — > oo leads to significant simpli- 
fication of the many body problem while retaining es- 
sential dynamical features of low dimensional situations. 
An important amount of studies has been done within 
this approach, mapping the lattice problem onto a self- 
consistent Anderson single impurity model. This map- 
ping becomes exact in the limit of infinite spatial di- 
mension. For a review, see Q. Most of this work 
has been done focusing on the original one-band Hub- 
bard model, which is the simplest model to describe the 
electron-electron interaction-driven metal-insulator tran- 
sition (MIT). 

However, orbital degeneracy in known to play a cru- 
cial role in correlated systems. The degeneracy of the d 
band is two in V2O3 [|| and three in LaTiOs. Orbital 
degree of freedom is relevant to explain some very in- 
teresting properties like colossal magneto- resistance 0, 
MIT in alkali-doped fuUerenes f^,0] , and for any physical 
property involving orbital ordering. 

Many theoretical approaches have been proposed to 
describe the effect of strong Coulomb interaction in 
systems with orbital degeneracy, using the slave-boson 
method p-[l2|, the variational method ||l^, and the 
limit of high spatial dimension. Some results of these 
works concerning the transition criteria and the order 
of the MIT are substantially different. An unified the- 
ory describing degenerate Hubbard model is still missing. 
Concerning high spatial dimension approaches, quan- 
tum Monte-Carlo (QMC) JT^ , p2l an d a generalized iter- 
ated perturbation theory (IPT) [|15[ have been proposed. 
Imaginary time results, high computational time and fun- 



damental difficulties at low temperature can be a limit for 
QMC calculations. The generalization of the IPT, devel- 
oped in Ref. iQ , presents the advantage of being able to 
deal with particle-hole asymmetric problems. However, 
this approach is essentially a perturbation theory with 
respect to the correlation strength, which is the highest 
energy scale of the system. The non-crossing approxima- 
tion (NCA) [|7| was shown to be in excellent agreement 
with quantum Monte Carlo calculations for the DMFT 
of the one-orbital case ||l8|,|9| , even far away from half- 
filling. 

Here, we present a new approach based on the NCA 
within the dynamical mean field theory (DMFT). 

The starting Hamiltonian is given by the two-orbital 
degenerate Hubbard model 
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where the sum < i,j > is the sum over nearest neigh- 
bor sites of a Bethe lattice and a,b — 1,2 is the band 
index, cf^^ (respectively Ciaa) denotes the creation (re- 
spectively annihilation) operator of an electron at the 
lattice site i with spin a and orbital index a and Uiaa is 
the occupation number per spin and per orbital. The on- 
site Coulomb repulsion U and the exchange parameter J 
are assumed to be independent of orbital. In addition 
we will neglect the last term in the Hamiltonian (|l|) and 
the hopping between different orbitals i°^ = —tSab- Do- 
ing such approximations leads to the same Hamiltonian 
studied previously by quantum Monte-Carlo (QMC) Q. 
Integrating out the fermionic degrees of freedom leads 
to the single-site effective action 
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Therefore, as in the one-band case, the large- -D version 
of the degenerate two-band Hubbard model is mapped 
onto an effective impurity model. The local correlated 
site is now orbitally degenerated, and two different ef- 
fective media have to be considered. The corresponding 
Hamiltonian reads 
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where the local part is 
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and the coupling with the effective medium is 
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The orbital dependent coupling with the effective 
medium is characterized by the hybridization 
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From the equation of motion of this effective Hamilto- 
nian, it follows that 

Ga.{uj)-^ - ^ - Erf - Saa(cc') - J%UJ) . 

Comparing this equation to the following property of 
Green's functions on a Bethe lattice 

we found the self-consistent set of equations 

In the following, we solve the impurity model of equa- 
tion (S) using the extended version of the non-crossing 
approximation presented in a previous work pO| . Propa- 
gators and self-energies of the sixteen local states \a, f3 > 
are introduced, where a (respectively /3) is the local oc- 
cupation of orbital a = 1 (respectively a = 2). Each 



local state |a, /? > is an eigenstate of the local part of the 
model Hamiltonian. The corresponding local model that 
have to be solved is schematically represented in figure H. 
The impurity site is coupled with two effective media by 
two effective dynamical hybridization J'^(w) and J^{uj). 
Local states |a,/3 > and |a',/3' > are coupled by NCA 
equations if one of the two following conditions are ful- 
filled :|noi(a) — noi{a')\ — 1 and f3 = (3' or a — a' and 
\no2{(3) — "o2(/3')| = 1- The first and the second cases 
involve respectively J^{uj) and j'^{u}). 
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FIG. 1. Two orbitals local model. The sixteen local states 
|q,/3 > are all coupled by the orbital dependent effective 
medium hybridization J'^ito) and J^luj). 

Here are the main results of our dynamical mean 
field theory. The NCA approach presents the advantage 
to provide directly real frequencies one-particle Green's 
functions. Local quantities like densities of states can 
therefore be computed for any regime of parameters. 
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FIG. 2. Density of states per spin and per orbital for 
[/ = 4 eV, i = 1 eV and r = 1000 K. Total occupation 
number is n = 0.987. The inset shows filling dependence of 
low energy excitations. Solid line (respectively dotted and 
dashed line) is the density of states around the Fermi level 
for n = 0.987 (respectively n = 0.992 and n = 0.995). 



NCA was in very good agreement with essentially ex- 
act quantum Monte-Carlo simulations for the one-band 
case. This approximation is then a good candidate for 
the investigation of the multi-band generalization of the 
Hubbard model. 

Useful and comprehensive information concerning elec- 
tronic dynamic processes are given by these densities of 
states, without any size effects due to small clusters and 
till very low temperature. 

Spin and orbital dependent densities of states paai^) 
are obtained from NCA Green's functions by : 
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Figure displays the density of states per spin and per 
orbital of the doubly degenerate Hubbard model for a to- 
tal occupation number n = 0.987. Both lower and upper 
Hubbard bands are present. Relative spectral weights be- 
tween the lower and the upper Hubbard band are given 
by the ratio 1/3. Quasiparticle excitations can be seen 
around the Fermi level. The origin of this low energy 
structure can be clarified in the context of NCA by ob- 
serving its temperature and filling dependence p9[ or by 
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FIG. 3. Comparison of the densities of states per spin and 
per orbital for [/ = 5 eV, i = 1 eV and T = 1000 K in the 
one-orbital case (dashed line) and the two-orbitals case (solid 
line) Total occupation number is n = 0.94 in both cases. The 
inset shows low energy excitations in both cases for the same 
parameters. 

reducing the Hilbert space of available local states onto 
the impurity pffl. 

The proximity of the n = \ metal insulator transition 
is confirmed by the disappearance of the low energy spec- 
tral weight. The inset shows the behavior of quasiparticle 
excitations just before the transition. Total occupation 
numbers are n — 0.987, n — 0.992 and n — 0.995. 

It is interesting to compare two-bands densities of 
states with one-band NCA results. This is shown in fig- 



ure ^, for [/ = 5 eV and ritot = 0.94. Spectral weights 
are in good agreement with previous approaches like the 
IPT one |13]. More precisely, we found that degeneracy 
has a strong influence only on high energy excitations : 
respective spectral weights of lower and upper Hubbard 
bands for ntot — 1 are 1 and 3 in the two-orbitals prob- 
lem, instead of 1 and 1 for the same filling in the one- 
orbital problem. However, NCA bands structures are 
strongly different to IPT ones. Fine structures displayed 
by IPT results are not present in our approach which is 
confirmed by recently obtained QMC spectral functions 
for the multi-orbital Hubbard model in infinite dimen- 
sion 11]. 

Figure Q shows the general evolution with respect to 
the total occupation number of electron density of states 
per spin and per orbital, for t/ = 5 eV, t = 1 eV and 
T — 1000 K. Each curve corresponds to a given Fermi 
level. Quasiparticle low energy excitations can be seen 
around the Fermi level. Total occupation numbers n are 
indicated for each density of states. For integer fillings 
n = 1, n = 2 and n = 3, the system is an insulator. It 
is metallic for any other filling. The interesting transfer 
of spectral weight from high to low (Fermi level) energy 
scale will be studied in details later. 
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FIG. 4. Densities of states per spin and per orbital for var- 
ious fillings. The correlation strength is [/ = 5 eV, bandwidth 
is t = 1 eV and temperature is T = 1000 K. For each curve, 
Fermi level is labeled by a black point. 

Note that effective mass of quasiparticles increases 
when approaching the Mott insulating state, which is 



consistent with experimental observations in filling con- 
trolled oxides family such as Sri_a;Laa;Ti03 [ pi| |. 

The Fermi level dependence of total occupancy is plot- 
ted in figure ^ These results are consistent with re- 
cent QMC simulations for the multi-orbital Hubbard 
model p2]. When the system in an insulator, as for 
U ~ A eV in the figure, a small displacement of the 
Fermi level in the insulating gap does not lead to a sig- 
nificant variation of the total occupation number. A plot 
of occupation number versus chemical potential displays 
a plateau as soon as a MIT is crossed. Our results show 
that three successive MIT at integer fillings occur in the 
system, n = 1 and n = 3 transitions are equivalent be- 
cause of the particle-hole symmetry. For a given U, the 
Mott gap of the n = 2 transition is smaller than the n = 1 
gap. This is in good agreement with previous calculations 
leading to a larger critical Uc for the n — 2 transition [Q . 
We found Uc{n = 2) « 3.5 and Uc{n = 1) « 3.1. 
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FIG. 5. Total occupation number with respect to Fermi 

level of d electrons. Parameters are U = 5 eV, f = 1 eV and 

temperature is T = 1000 K. The inset shows the occupation 

up to the transition value n = 1, for [/ = 4 eV. 

In summary, we have proposed a new approach based 
on the dynamical mean field theory for the doubly de- 
generate Hubbard model, combined with a generalization 
of the non-crossing approximation in order to solve the 
equivalent self-consistent local impurity problem. This 
impurity problem consists in a doubly degenerate local 
correlated site (sixteen local states) embedded in a con- 
duction band which is orbital dependent. 

We obtained a description of the integer filling metal 
insulator transitions consistent with previous theoretical 
works using slave-boson formalism, variational method 
and quantum Monte-Carlo dynamical mean field theory. 
One of the advantages of our approach is to provide the 
full spectrum of excitations for any regime of parame- 
ters and a clear interpretation of the quasiparticle na- 



ture. Besides, many extensions of the model Hamiltonian 
can be done within this framework. In future calcula- 
tions, we will increase the degeneracy (some oxides, like 
Sri_a;Laa;Ti03 are triply degenerate) and we will take 
explicitly into account the exchange interaction by intro- 
ducing a non zero Hund coupling parameter J. We will 
also extend our analysis to thermodynamic and transport 
properties. 
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